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^STRACT 


In this thesis an attempt is made to model the signal 
received from a target, such as a submarine, by an active 
SODaR, such as a homing torpedo. For an iso-speed shallow 
water channel, the channel exhibits discrete multipath 
scattering. For typical LODAR beam patterns and normally 
encountered boundary reflection losses, the number of paths 
with significant energies contributing to the received 
signal is small. 

A characteristic feature of the above channel is seen 
to be the common doppler shift in the composite received 
signal due to target motion. This doppler shift is taken 
to be the average of all the doppler shifts of contributing 
rays as doppler ratio differences among paths are seen to 
be small. 

Modelling the attenuation coefficients as Zeromean 
complex Gaussian random variables and modelling the random 
travel times with appropriate density functions, the 
optimum receiver is found to be too complex to be of practical 
interest. If perfect knowledge of arrival times is assumed, 
the optimum receiver is of simple nature. A conventional 
receiver under this characterization suffers no marked degra- 
dation. With an appropriate choice of envelope of the trans- 
mitting pulse, the optimum receiver implementation is simpler. 
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CHAPTER 1 
INTRODUCTION 

The normally used receiver for the detection of a 
target in presence of white noise by an active SODAR (S und 
Detection And Ranging) is a correlation or a matched filter 
receiver. For such receivers the channel is considered to 
be dispersionless and unbounded, and the target is modelled 
as a slowly fluctuating point target. Then the return signal 
in an active SODAR is a replica of the transmitted pulse 
appropriately weighted and time delayed. Even if coloured 
noise is present, the same techniques used above suffice and 
it is generally observed that a suitable choice of trans- 
mitting waveforms will improve the receiver performance, 

? For bounded channels, there was no special mention 

on the performance changes for the above receivers. In this 
thesis we attempt to restructure the target and channel 
models when the channel is bounded and suggest some suitable 
receiver structures, 

A bounded channel is assumed to have well defined 
closely spaced boundaries, (an example could be the shallow 
water channel where the source is operating at high frequencies) 
and such boundaries are said to produce specular reflection. 
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at high frequencies it is appropriate to characterize the 
reflection process to be regular and deterministic. Now 
in such a channel v/hen the source/receiver and target are at 
fixed locations, the acoustic transmission may bo described 
by ray theory and the channel is said to have a discrete 
multipath structure. Hence the received signal could be 
a super-position of v\?eighted and time delayed replices of 
the transmitted pulse arriving at receiver after the 
'reflection’ at the target. In addition, a moving target 
will also introduce doppler shift on carrier frequency along 
each path. Hence the net effect of the channel and a target, 
which is modelled to have zeromean complex Gaussian r ef lection 
coefficient, is the spread of received energy over ambiguity 
plane. The a aount of spreading of the received signal over 
ambiguity plane may be found by using ray tracing programmes 
even for a channel in which the sound transmission is a 
function of depth [ij. 

It is observed that the general interest in a multipath 
environment in a bounded channel is to investigate the 
coherence degradation between signals received at two 
remotely located sensors [2,3], This scattering phenomena 
is some times applied in passive tracking of a moving 
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Source in shallo’-v waters [4,5]. An elegant model for the 
channel is suggested by Jacyna, Jacobson and Clark [5], 
where a passive source is tracked by two sensors. 

In this thesis vje adopt the Jacyna, Jacobson and 
Clark's model for an active SCDAR to detect a moving 
target. 

In certain occasions it is observed that special 
signal v/aveforms like coded pulse sequences or analog 
modulated v;aveforms (LFM etc.) are used in channels where 
thermal noise is a predominent component, for detection 
of targets. In this thesis we attempt to analyse the 
requirements of such waveforms when the channel is bounded. 

In Chapter 1 of this thesis the problem under study 
is proposed. In Chapter 2 , we first review various concepts 
that are involved for computing the parameters needed in a 
multipath received signal formulation, and adopt a model for 
detecting a moving target in a shallow water channel for an 
active oODAR. In Chapter 3, we attempt to study the optimum 
receiver for the received signal characterized by random 
attenuation coefficients and random time of arrivals. V»e 
then observe the simplicity in optimum receiver design 
having a knowledge of time of arrivals. For this situation 
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we consider a conventional receiver and study the amount 
of its degradation coraparod to the optimum receiver. The 
effect of choice of transmitting waveforms are studied. 

In the end we discuss the results. In Chapter 4, v-ie 
conclude the discussion v\/ith some suggestions for future 
work. 



CHAPTER 2 


IviODELLINS OF S'rinLLOV/ WATER TARGETS 

2,1 INTRODUCTION: 

In this chapter we are interested in developing 
a model for an under Vi/ater SODAR target. Our concern is 
to characterize the signal received by an under v/ater weapon 
system when it transmits a signal that is ’reflected' from 
a moving target. Such a situation for example, is encoun- 
tered when an active torpedo homes onto a submarine moving 
in continental shelfs characterized by shallow waters. 

The received signal features are dependent on -the 
channel parameters which modify the acoustic propagation of 
the transmitted signal. Boundary conditions and channel 
parameters may have complex influence on the acoustic 
propagation. In Sec. 2.2, we consider the effect of' various 
parameteis the medium on wave propagation and the 
relevance of ray acoustics to the present problem. In Sec* 
2,3, we give expressions for travel time, spreading loss and 
reflection associated with an acoustic wave. In Sec. 2,4, 
we consider the multipath structure of the channel. The 
influence of a moving target and a general treatment to 
compute the return signal in a moving target environment is 
studied in Sec^, 2,5. Apd finally in-Sec, 2*6 we derive^an 
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expression for the return signal for the above specified 
example and discuss the results obtained. Our account of 
the target model, in particular the influence of the source 
motion on the received acoustic field, closely follows the 
work of G.ivi, Jacyna, M.J. Jacobson and J.Q. Clark [5], 

2.2 U^JDER- water ACOUSTIC PROPAGATION PLINCIPL&Ss 

The acoustic wave propagation in an underv\/ater 
channel is influenced by regular and random Inhomogeneities 
of the medium [ d]-.- Refraction, duct propagation etc. are 
the phenomena due to regular inhomogeneities. The regular 
inhomogeneities refer to the spatial and temporal variations 
in the mean characteristics like temperature, salinity and 
pressure etc, in the medium. In ocean media these effects 
are considerable in the direction of depth and the ocean 
medium is often said to have a layered structure. The 
random inhomogeneities on the other hand, cause phenomena 
like scattering or reverbaration, Which of the above two 
inhomogeneities is dominant depends on the specific situation 
under study. In this thesis we .are mainly interested in 
regular inhomogeneities only. 

The regular variations in temperature, salinity and 
pressure influence the magnitude of the velocity of propa- 
gation of the acoustic wave and consequently the path 
traversed by the acoustic wave. 
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ye now briefly consider the basic principles which 
govern the acoustic propagation in the medium* 

The propagation of sound in an under v/ater channel 
can be described mathematically by the solution of the wave 
equation using the appropriate boundary and initial condi- 
tions for a particular problem. In a homogeneous and 
quiescent medium the velocity potential 0 at a point on the 
v;ave obeys the wave equation given by [ 7 ] 

0 - = 0 (2.1) 

c dt 

where *c’ is the velocity of propagation of the v/ave and 
c is given as 
c = 

where K is bulk modulus and is density of the medium. 

The velocity of the particle on the wave is given by 

V = 

and the pressure developed due to this velocity is given 
by 



where is the density of the medium. 
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Under the assumption of small perturbations in an ambient 
state, the above wave equation can easily be seen to follow 
from? 

th e principle of c o nserv at ion of mas st This principle 
states that in a fluid the net rate of increase of mass per 
unit volume is equal to the rate of change of density. If 
P(x,t) is the density of fluid at a point X and time t, 
v(x,t) is the fluid velocity at that'upoint, the equation for 
conservation of mass is easily shown to be given by 

If -K V' iP^) - 0 

ii) The Euler’s equation of mot ion fox a fluid : This equation 
states that the mass multiplied by acceleration of a fluid 
particle is equal to the net apparent forces exerted on it 
by it*s environment and by external bodies. Neglecting the 

i 

body forces such as gravity per unit mass, the equation of 
motion is written as 

P ■ ^ = -Vp 

where p(x,t) is force per unit area or pressure. 

The equation of sta te; It describes the relation 
between the pressure of a fluid, and its density and temperature. 
If an excess pressure developed due to an acoustic wave is 
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p, then this excess pressure is proportional to the fra- 
ctional change in density (Ap)* Hence v;e write 
Ap = K-A/°w/here K is the bulk modulus of the medium. 

There are two theoretical approaches to the solution 
of wave equation (2.1). One is called the normal mode 
theory [8 in which the wave propagation is described 
in terms of characteristic functions called normal modes, 
each of v/hich is a solution of equation (2.1). In a bounded 
medium the standing wave pattern occurs due to super position 
of waves travelling in opposite directions towards boundaries. 
If H is the depth of the channel, then for all v^avelength 
equal to where j = 1,2,3,..., a solution to wave 

equation is said to exist. In ocean media, the standing wave 
pattern occurs along the z— axis (depth) and the wave propa- 
gates along the x-axis (horizontal axis). The normal modes 
are combined additively to satisfy the boundary and initial 
conditions. With normal mode theory we obtain at a point 
and at a given time the combined effect of all wave paths. 
This analysis is relatively easy to handle when the number 
of modes is small. If the point at which the effect of the 
acoustic source is computed, is at a large distance (of the 
order of 50 I<m) or the source frequency is small (less than 
1000 Hz) the resulting number of modes will be small and 
normal mode theory may be applied. 
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The other approach to solving equation (2.1) is by 
ray theory [9 ]. Ray theory represents the acoustic field 
at a point as a sum of ray contributions, each ray emanating 
from the source and reaching the point after a certain 
number of boundary reflections. In media in which the sound 
velocity changes appreciably over a small distance (say in 
the order a wavelength of the sound v;ave), ray theory is not 
applicable. Hence ray theory may be seen to be applicable 
to high frequencies or short wavelengths in a practical 
situation, liay tracing can be done using suitable computer 
programes for the energy loss, time of travel along each 
ray path [ 10 ], 

An underwater weapon generally has the following 
operating specif icationst 

i) Source frequency range : 10 - 100 KHz 

ii) Acquisition range : 2-3 KIvi 

iii) Operating channel depths : 200 - 400 meters. 

These factors allow us to use ray theory to analyse 
the received signal characteristics of an underwater 
weapon, 

2 . 3 RAY ACOUSTICS: 

We consider the wave equation in two dimensions, say 
X and z where x signifies the horizontal axis and z signifies 
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the depth, in a ocean medium. The positive direction of z 
is tovi?ards increasing depth. Let the acoustic source be 
producing a harmonic disturbance, :;e consider a homogeneous 
medium in which velocity of propagation of the acoustic 
wave 'c* is constant. Then the wave equation (2,1) is 
written as 



consider 

0(x,z,t) =s'Y(z) exp j(+ax~<.ot) (2.3) 

to be a solution of equation (2,2), 

hquation (2,3) describes plane v/aves propagating in the 
direction of increasing x (+ sign) or decreasing x(- sign). 
The phase velocity of the wave is given by v and is 
equal to ~ , 


Substituting equation (2.3) in (2,2) gives the differential 
equation 



2 

yy 


0 


where 



a 


2 ,v 
a (“j 

c 


1 ) 


(2*4) 


(2.5) 
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Consider a solution of (2.4) as 

(z) = (2.6) 

where A and 3 are constants of integration, 

Coiiibining (2,6) and (2.3), we have 

«l(x,z,t) = A + B (2.7) 

where the first term describes a wave travelling in the 

direction of increasing z, i,e,, propagating down-wards 

and second term corresponds to a wave propagating upwards, 

2 

we write equation (2,5) as =s and define 

c 

4 ( 2 . 8 ) 

c 

where k is the wave number in the direction of propagation. 

It is the magnitude of the vector K whose direction is normal 
to the wave front i.e., K defines the ray direction. 

If we consider a plane wave described by an 
equation 

0(x,z,t) = .. 

then this plane wave is moving in a direction K whose magnitude 
k is given by or inturn we can say that the plane wave 

is having a and y as its x and z components. 

Consider Fig. 2,1 and if Q is the angle the direction 
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of plane wave makes with z-axis, then the horizontal component 
a = k sin0 anct vertical component y = cosQ. Fox a homo- 
geneous medium we can see K is constant. That is the wave 
will , not alter it’s direction and it’s magnitude remains 
constant as it propagates. 

Now v;e derive a solution to the wave equation (2.2) in 
a medium where c is not a constant. 


For a source producing harmonic disturbance consider 
a solution to the wave equation (2,2) as 

|i!(x,z,t) =A(x,z) (2.9) 

where b, a function of x and z to be determined, we substitute 
(2,9) into wave equation (2,2) and separate the real and the 
imaginary parts. Ignoring the time dependence, we have 


A + K^A = O 

a nd 


2 V A .V ^ A = 0 

Rewriting (2.10) and (2.11) as 

2 , 


+ C||)^ - k2 




and 


V^S+I VA.'eS = 0 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 
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then equation (2.12) implies 


f oS\2 


/dS^2 

(■gy) 




(2.14) 


Hquation (2.14) is called Eikonal equation (corresponding to 
the wave equation) and S(x,z) must satisfy (2,14). S(x,z), 

O 

leading to a solution to the wave equation, if -■ VA.*^ S is 

A 

Small in equation (2.13) must also satisfy 


= 0 (2.15) 

The general form of the solution of equations (2.14) and 
(2.15) will be 


S(x,y) = + ir(K^ ~ ^1^)* ^ (2.16) 

where c^^ is a constant to be determined. 

Equation (2.16) describes a plane wave front in which K is 
the wave number. It is the magnitude of a vector K whose 
direction describes the ray path. Let K make an angle 0 
with horizontal axis (x-axis), then we can show that 
Cj^ = K cosa. This implies that c^ is the direction cosine 
of K along X'-axis, 

Having obtained S(x,z), we substitute equation (2,9) 
now into wave equation and obtain a second order differential 
equation in A(x,z), Solving for A(x,z), we obtain 0(x,z) 
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Conditions 


.,2 


^ — "SiO and ^ A .^3 is small are satisfied, if 
the fractional changes in medium parameters like velocity are 
small over a wavelength. These conditions are generally met 
when the source operates at high frequencies. 


The solution, (2,16) of eikonal equation is seen to be 
a plane wave with Kk ~ , If velocity of sound ’c‘ is a function 
of X and z, the gradient of the wave front and therefore of the 
ray path change with x and z. Now we derive a condition which 
the ray path has to satisfy at each point. 


Each point on the ray satisfies the eikonal equation 

(||)2 + (||)2 = k 2 = 4 • 

c 

An infinitesimal ray length can be described as 


dd h. kdl = k(dx + dz) 


1/2 


(2.17) 


using Fermet*s principle (which states that rays in physical 

space are paths of stationary time i.e., for small variations 

in path length between two fixed points, the resulting varia- 

^2 

tions on the travel time integral / dt is zero); if the 

variation introduced on the path length between two fixed points 

vw have ^ f do = 0, then we have from equation (2.17), for a)»0 

S / K df =<?/ O (2.iS) 

If we consider the velocity of wave ’c' to be a function 

of z only, which is true generally in ocean media, then 
dx 

writing ^2 ~ Fz (2.18) implies 

2 

ctry 


V". X 

Sf 


+ 1 ^ 


dz = 0 


(2.19) 
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Applying Eulex-Lagrangss formula, to equation (2.19) 


we have 


d _ . dF 

f 


dir 57 " = ° 


( 2 . 20 ) 


where F 


“ETTT- 


dF 


Equation (2.20) implies - = constant and w'e have 


1 

c 


z _ constant (=a) (2.21) 


r +1} 


From Figure 2.2, if 0 is the angle the ray makes with x-axis 
then 


X = -^1 = cot0 
Z GZ 


(2,22) 


Therefore (2. 21) implies, cos0/c( z) = a 

COS0 _ , 

cIYJ - ® 


(2.23) 


Equation (2.23) is called Snell's law and it is to be 
satisfied at each point on the ray. 


The eikonal equation in 3-u can be written as 
5yJ 


(i)^ + (f)^ + (il)^ = k2 


If we assume the velocity of the wave *c' is a function of 
one co-ordinate only, namely z-coordinate, then for a ray 
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path, the direction cosines dS/dx and dS/dy along x and y 
co-ordinates must be constant at all points. This implies 

= constant and hence ray path is confined to a plane 
normal to xy plane, oo with no loss of generality we 
assume the ray plane to be xz plane. 


2,2,2 Computation of Travel Time and bpreading Loss: 


',Ye have seen in a medium in which 'c' is proportional 
to only z-co-ordinate, the rays are to be confined to xz 
plane only, b'e have also seen that at each point on the ray 
the Snell* s law is satisfied, he now compute the horizontal 
distance (along x-axis) a ray traverses when it travels from 
a given depth (say z = O) to another depth (say z = z). 
Equation (2*22) implies dx = cot© dz, and the horizontal 
distance traversed is given as 

z 

x = / cot © dz (2.24) 

o 





dz 


(from equation (2,23)) 


If *0* is a linear function of zj i,e. writing c = pz, where 
p is the gradient, then equation (2,24) is evaluated and 
we have 

X = const, -h (1-a^p^z^)^/^ (2.25) 
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Define HA and ^ const,, the (2,25) implies 

(X - X^)^ + (2.26) 

Hquation (2,26) implies that the acoustic ray traverses 

a circular path if *c' is a linear function of depth. The 

radius of the circular path is and the origin will lie 

ap 

at (Xj^,0) where z = 0 implies the depth where 'c' is 
extrapolated to zero. Fig. 2.3 illustrates the above 
statement. 


Let 0 be the angle of inclination the ray makes at 

the source with the horizontal axis. Let c be the velocity 

o 

of propagation at the source. From equation (2,23) we have 


cos 0 

“cTzT 


a 


1 c 

Radius of curvature R = “~ =s 

ap p COS0 

and the origin lies at a depth v^here velocity profile 

c 

extrapolates to zero or it is at a depth from the location 

of the source. Therefore centre of the circle is at (X,,0) 
c 

where X. = —3— tan©, 

1 P 

'v/e define travel time along the path is equal to 


r’S dz 

-Z- where dS = as may seen from Fig. 2.2, 

Therefore if c is a linear function of z (i.e. , c = pz) then 
travel time ’t* is given as 
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t = 


p2f(l-a2p^z^J) 


( 2 . 27 ) 


Transmission loss alonc; a ray path is a sum of 
spreading loss and attenuation loss. Attenuation loss 
includes losses -due to absorption and scattering. In this 
discussion vve ignore attenuation losses. Spreading loss is 
a geometrical effect representing the regular weakening of 
signal as it spreads outward from the source. In shallow 
water channels the medium is bounded by two surfaces and 
the spreading loss is assumed to be cylindrical loss and 
is given as [> ] 



where 

p is the source energy per unit time per unit solid 
angle , 

dJ-K is the area swept oaf by the wave surface, normal 
to the ray at a horizontal distance then 



X sinQ dx i~l 
p COS0^ 00^ ' 


( 2 . 28 ) 


whore 0^ is the angle the ray path makes with the horizontal 
at the source and 9 corresponds to the angle the ray makes with 
the horizontal at the point of interest. 
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2«3*2 Reflection Coefficient: 

At the interface separating two media, the ray has 
to satisfy the following ooundary conditions: 

i) the vertical particle velocity must be continuous. 

d0 40, 

-a-z- = -3f 


ii) continuity of pressure across the boundary is to be 
maintained, i*e. 


p P 

' 1 


d02 

Tt = * 2 ■"St 


•(2.30) 


where P ^ and densities of the media. 


Define 0^^ = A exp (ax + yz - wt) as incident ray and 

0 = B exp (ax - yz - ojt) as the reflected ray. Refer 

ye 

Fig, 2.4, 


Substituting above equations into boundary conditions, 
the reflection coefficient turns out to be 


B 

A 


f^2 

(^1 




(c-l^/c^)^- cos'^e^'j 

^ I 

1 - cos 9, -* 


Vv- 


Pl 




(c,/c„)^ - COS^©T I Vl- 


'r 12 

1 - cos 9, 


f 


(2.31) 
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where and axe the velocity of propagation and 
density of ocean medium and c^ and corresponds to same 
quantities in the second medium. 


Consider the second medium to be air and the interface 
to be ocean surface. No»v the interface separates two media 
having large difference in their acoustic impedance. Acoustic 
impedance is defined as P c, and we have c^ = 1,5x10^ g/ 

(cm) sec, and ^2^2 ~ gm/(cm)^sec. (for air), 

So equation (2.3l) implies that the reflection coefficient is 
unity and the reflected ray will have it’s phase advanced 
by 180°. 

On the other hand at the interface separating the 
sea and the sea bottom with c^ as the parameters corres- 
ponding to the sea bottom, v;e have Cj^ < c^ and ^ p2* 

ive define critical angle 9^ as 

c 

c, 

cos (2,32) 


Now from equation (2,3l) and (2,32) we see that the 

reflection coefficient v/ill increase from a value equal 
^ c ^ c 

to for © = 90° to a value, unity at angle of 

1^1 

inclination © = ©_. A typical curve for reflectivity vs. 
inclination angle is shown in Fig. 2.5, For the angle of 
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inclination © from 0 to the reflection coefficient is 
purely imaginary with magnitude of the reflection coefficient 
being unity and the phase change is determined by inclination 
angle. The phase angle is given as 


6=2 tan 


-1 


fcos^e-Cc^/c^)^ 

T2 


'/o 


rcos ©-.( C ,/ C ^) 1 

/ ^,..1 / ( 2 . 33 ) 

L 1 - cos^e J 


and ^ is 0° at the critical angle and increases to 180° 
at 0 = 0°. 

2.3,3 Ray Displacement; 

In general, the fields computed by the ray theory 
differ from those computed by the normal mode theory. Tindle 
and Bold [l2j account for this discrepancy by the beam 
displacement suffered by rays on total reflection from the 
bottom. Total reflection means that for rays v/hose inclin- 
ation angle are less than critical angle, the reflection 
coefficient amplitude is unity and has only phase change. 

And this additional phase is considered as extra path. When 
the extra path lengths and travel times are incorporated into 
the computation of acoustic fields, the two approaches were 
seen to give very nearly the same result, 

Brokhovskikh [13] has shownm that the lateral displa- 
cement A of a bounded beam for above situation is given 

2i= - de/da 
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where €■ is the phase of the reflection coefficient given in 
(2.33) and is 

6. = 2 tan 

and a is the horizontal component of the wave number and 

is given as cc = cos 0, where w is the operating 

^1 

frequency of the source. 

The time displacement 'X associated with beam displacement 
is given as 

“f — 

“ "'STT 


L 


9 , / \ 2 

cos ©-(Cj^/c^) "“j 

^ ^ • ..rr-- j 

l—cos © 


ivaluating and T , v/e will have 

2 . 2 


^ ^1*^2 ® 


^2 2 ^ 

Pi <=2 > 


and 

r = 


Pi 

rcT^^^cos^ - (Fi^c^^cos\ H-P^^c^^sin^e-F^^c^^) 


■or evaluating A and X we make use of direction cosines 


given as 


w 


Y, = sin© 
^ '^1 


w / i, \l/2 

Y 2 = c~ ^ 


a = — - cos© 
*=1 


0 being the angle the ray makes with horizontal. 



e 
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Thus is a function of frequency. The value of 

A is much less than a meter per reflection for frequencies 
that are of interest to, this thesis. For moderately long 
ranges, the additional path to be considered is negligibly 
small and so the ray displacement contributions are ignored 
in further discussions. For example 
for f = 30 KHz 
we have A = 5 cm* 

2,4 fv'iULTIPATH STRUCTURE OF THE Ca-JTHSL: 

Up to this point in this section, we developed various 
equations to compute ray path parameters. In this section v/e 
further investigate the nature of the ray paths. 

In Section 2,3 we have seen that for a medium in 
v^hich velocity profile is a complex function of depth, every 
point on the ray has to satisfy the Snell's law and a ray path 
is confined to a plane noraal to xy plane. Thus the various 
ray paths can be traced in principle. Depending on the 
complexity of the profile and the permissible resolution 
between paths, for each ray starting from the source, vie 
consider the range over which the velocity 'c' can be 
considered to be a constant or a linear function of depth. 
Accordingly v/e approximate the ray path either with a straight 
line or with an arc of circle, trace the path considering 
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the successive segments over which velocity profile is 
approximated. Thus we can draw all the ray paths that 
eminate from the source. Depending on the location of target, 
we can compute the number of rays that intercept a target. 

For a given configuration of source and target locations and 
channel depth, wc may be able to plot in finitely many 
paths connecting the source and the target. 

Fox ease in handling, the following nomenclature is 
often adopted in classifying various rays. Among the group 
of rays that are obtained, we have a direct ray which is the 
line of sight ray. The remaining rays are ordered according 
to the number (n) of bottom reflections the ray has suffered. 

This number is called the order of the ray. Now n can have 
all integer values starting from zero. For n=0, i.e. for 
the zeroth order ray we have only one surface reflection. For 
each order of bottom reflection Vi/ith n > 1, we can associate 
a set of four paths; one of them will be having n+1 surface 
reflections, another one will have n— 1 surface reflection and 
two will have n surface reflections. The direct ray, and the 
rays having order n=0,l,2 are illustrated in Fig. 2.6. 

we assume the velocity of sound c at a point in the 
medium is linearly proportional to the depth of the point in 
the medium. Depending on the positivity or negativity of the 
gradient, location of the source and the initial ray inclination. 
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each ray travels in one of the four characteristic v;ays. 

These categories axe: 

i) the ray which gets reflected at both the interfaces 

ii) the ray vjhich gets refracted at both the interfaces, i,e, 
the ray will not reach the interfaces but curves away 
due to the linear dependance of c on depth. 

iii) the ray v^/hich gets reflected at surface and gets 
refracted at bottom (RSk) 

iv) the ray which gets refracted at surface and gets 
reflected at bottom (RBR,)» 

These four types of rays are illustrated in the 
Figure 2.7. 

M.J. Jacobson [14] considered the ray paths under 
RSiT and RBI. categories and gave formulae to find the number 
of paths arriving at a given point from a source, path 
losses, path travel time and the initial inclination angles. 

We now consider the approximations v/e make and the 
simplifications that result on ray path tracing in the channel 
of our interest. 

We assume an iso-velocity shallow water channel 
which is having flat boundaries. Therefore ray paths 
travel in straight lines irrespective of the location of 
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the Source and all rays belong to the category (i). The 
rays will have specular reflection at the boundaries. The 
direct ray will have spherical spreading and due to good 
reflectivity conditions assumed for boundaries, rest of the 
rays will have cylindrical spreading. Hence received energy 
decreases in proportion to the square of distance for 
targets at closer ranges and decreases in proportion to the 
distance for longer ranges. The under water weapons v;ili be 
generally active in the ranges ranging from 500 meters to 
2500 meters. Therefore on the average vje assume the loss is 
due to cylindrical spreading. 

From the discrete ray structure shown in Fig, 2,6, we 
can see that the direct ray arrives first at the target and 
is followed by either S ray ox n=l ray depending on the 
location of source and target. Then the higher order rays 
will arrive. As the order of the ray 'n* increases, the ray 
inclination at the source will become steeper. If the 
initial inclination is steeper than the critical angle (©^) 
defined in sec. 2.3,2, theray suffers reflection loss in 
magnitude. The number of reys that are to be considered are 
decided by the receiver and transmitter beam patterns. Among 
them we may neglect the rays that have inclination angles 
steeper tl^an critical angle and undergo more bottom reflections. 
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On the average for a typical values of ^2 ® 

ray suffers approximately 10 dB energy loss per bottom 
reflection of inclination angle exceeds critical angle. 

Till nov/ we consider the boundaries to be flat, but 
in a physical situation the bottom structure may have either 
positive slope or negative slope. The ray tracing can be 
carried out with a simple modification of the above procedure. 
This is discussed by G.V. .-.nand and P, Balasubramanian in 
their report [l5j, 

2,5 EFFECT OF TARGET MOTION ON RECEIVED SIGNAL CHARACTERISTICS: 

Using ray theory, we have seen that shallov; water 
channel gives rise to a discrete multipath structure to the 
received signal. The received signal energy now is obtained 
by superposition of the energies arriving along these paths 
after getting ' ref lect^jd’ at the target. If ‘n’ significant 
paths are to be considered each of the path energies may be 
characterized by an attenuation coefficient (A^), arrival 
time (T^) and initial path inclination angle (©j^) . 

We now consider the effect of target motion, A 
moving target introduces a doppler shift in the frequency of 
the received signal with respect to the transmitted signal 
frequency. If V is the velocity associated with the target. 
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the doppier shift is defined in terms of the component of 
the velocity along the direct ray, he assume direct ray 
path length is large compared to the target movement over 
the time interval of interest and hence the only effect 
on the received signal is changing of phase with respect 
to receiver time frame ’t’ or equivalently doppler shift in 
the transmitted frequency. 


he may now define a doppler shift (k^^) along each 
path by considering the velocity component of the target 
along that path and it may be given as 


k 


n 


= 2 




where 9 ^ is the velocity component of target motion 
along nth ray. 

We now argue that difference between the doppler 
shifts of any two discrete paths is small as we view that 
the doppler shifts are due to the changes in the path 
lengths due to target motion. Hence v/e propose a common 
doppler shift for the total received signal which is inde- 
pendent of the order of the discrete ray path, vie do this 
first by computing the phase of the received signal along 
each path and then showing that the change in phase of the 
combined signal is independent of the order of the constituent 
rays by using some valid approximations. 
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We first write an expression for the combined 
received signal in terms of it’s constituent paths and 
for a general target motion. 

Consider the geometry of the system as shown in 
the Fig, 2,8, 

Let the transreceiver (TR) be located at origin ’O’, 
Suppose that the target at any instant of time follows a 
path a with arbitrary velocity ’c', sound velocity could 
be a function of time. Target depth may change v/ith time. 
Let the instantaneous position of the target (T ) be given 
by P, Then OP is the slant range and is denoted as R(t), 
Consider the plane X’Z \>/hich contains the directed line OP. 
The projection of OP on OX’ is OP’ and is called horizontal 
range. The target is assumed to follow a path cr as shov/n 
and V(t) = is the velocity vector, IVe assumed that 

’n’ discrete ray paths erninate from transmitter, get 
’reflected’ at the target and arrive back at the receiver, 
vie also assume that the return path from target traces 
the corresponding forward path. Let 0^ be the angle nth 
ray (PRj^) which leaves the target, makes with OX' , The 
convention for the sign of angle we adopt isj 0^^ is 
positive if the ray travels in positive direction of 
z-axis (i.e,, towards increasing depth) and 0^ is negative if 
the ray travels towards decreasing z. 
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l'i8 have seen in Sec. 2.3.4 that the first ray to arrive 
at the receiver after reflection at target is direct ray. 

It is followed by rays having bottom reflections of order 'n' 
zero and onwards. Since the target is moving, the instanta- 
neous position of P follows the path and accordingly the 
vertical planes containing different rays may be different. 

Since the target is moving, the signal associated 
with nth ray at the receiver time ’t' corresponds to a pulse 
reflected at a target position at an earlier time. For nth 
ray, let the ray leave the transmitter at time ’t^’ and reach 
the receiver at time ’t*. Then we have 

tn = t - 2 T^(t) (2.34) 

where is one-way travel time. 

Here we see t to be direct function of t and so we may 
n ' 

write 

= F^(t) (2.35) 

Let the transmitteo signal is rectangular pulse with 
carrier frequency 
write, 

s(t) = Pve[Hxp (jm^t)j 0 < t £ T (2,36) 

Now the received signal along nth ray is written as 
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3rnCt}=Re[A^(t-Tn(t))exp jCco^(t-2 t)+S^( t) ) ] 

(2.37) 

where A^(.) is the attenuation coefficient along the path. 
This includes path losses and reflection coefficient of the 
target "the phase introduced due to boundary 

reflections and also includes the phase of complex 

reflection coefficientof the target. The nth ray may have 
N bottom reflections and m surface reflections along it* s 
path. Then 

S^(t) = 't n t) (2*38) 

v/here ^ is the phase change per one bottom reflection* Each 
surface reflection introduces a phase change of n. In 
principle number of reflections and phase change due to 
bottcam reflection could be functions of time, 

A^(,) is given as 

\(t) = ia„| L~^(t) [Bj^(t)]'^'^^^^ (2.39) 

where L~^(t) is two way spreading loss and we assume the 
path length difference in forward and backward paths is 
negligible over the analysis interval, 

Bj^(t) is the path loss due to one bottom 
reflection ia^j is the modulus value of complex reflection 
co-efficient. 
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Now we want to represent the received signal as a 
complex envelope with an associated carrier frequency 
First v/e write 


r(t) = Re[£ /^(t) exp j | o)^t 


= Re[E A^(t) exp j (a)^t-0^{t)] 

(2.40) 

where 0^(t) =* -2o>^ T^(t) + S^{t) 

(2.41) 

Then we write 


i(t) = A(t) exp j(o)^t+0(t) 

(2.42) 

Therefore, 


A(t) = U Sin0^(t)]^ + u A^(t) 

cos 0^(t)]2 


(2.43) 

and 


, Z A^(t) sin 0_(t) 

0(t) - tan — - 

Z A (t) cos 0 (t) 

(2.44) 


We now introduce some assumptions and approximations 
and derive the composite signal for the case of interest 
to us. 

2.5 iSO-SPtED CHANNEL* 

We assume that the velocity of propagation ’c’ is a 
constant throughout the medium. The channel is bounded by 
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horizontal boundaries. The target is assumed to follovj a 
horizontal straight line path with constant speed v. 

As shown in Fig, 2.9, let P be the location of the 
transreceiver and Q be the location of target at t=0. how 
Pv>{ gives the slant range at t=0» Let PC<' Rt.^q de 
the horizontal range at t=0, and y be the angle which the 
target path projected on to the horizontal plane makes with 
pc,’ at t=50. ri is the instantaneous position of target and 
R(t) is the instantaneous slant range. It’s projection on 
XY plane is Rj^(t). 

From Fig, 2.10 we observe, transreceiver ‘P’ is at 

a depth of H„, target is at a depth of H-j. and the channel 
ri * 

depth is H. Vertical height of PQ is *n'. As we have 
seen previously, each order of reflection is associated with 
a set of four rays. Each ray is characterized by its initial 
inclination. If 9;^^ is the angle mode by the nth ray, we 
follow a convention and say ^ a = -fl, if the ray 

leaves the source downward; else < 0 and a = -1, and 
also p = +■! if the ray reaches target/receiver from above, 
else p = ~1. This convention helps in calculating path 
lengths easily. Some typical ray oaths in a moving target 
environment are shown in Fig, 2*10. Here the rays may not 


lie in the same plane 
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First we calculate the ray path length (one way) for 
N bottom reflections assuming target is stationary. For 
this we can use the method of images of the source that 
contribute energy at the target. Each ray may be identified 
by considering the number of times it crosses the surface 
and bottom or their images. From the Fig. 2.11 we see that 

= nth ray path length 

= + p (2.45) 

where N is the number of bottom reflections and Fi is the 

o 

direct ray length. 

In Fig, 2,11, the line segment from pto 0 indicates 
the 0th order ray path (surface reflection only) and line 
segments from pto denote the set of 1st 

order ray paths, Fiay path lengths can be easily computed 
from the geometry. For example, Ray corresponds a path 
having two surface reflections and one bottom reflection. Its 
length is given as £r^^+(2H + H ^ * 

Therefore equation (2.45) gives nth ray path length 

/ 

as for a target at rest. Then one way travel time 

is given as 



(2.46) 
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Fxoin Fig, 2,9 we can compute the path length fox nth xay 
fox a moving taxget. 


= CC + * 2vt„ Rho cosr]l/2 


= R„Cl- JV cosy]^'^^ 


no 

where h = H 

T o 


n o' 


il 


(2*47) 


From equation (2»34) we write now 




n 


= t-2c'^i{v RqV t^(l- ^ 


cosy 


+ c2 T 2 i 1/2 
no •' 


= t— 2 [M^tj^^+2ri.^ c ^ M cosy 

*^0 


4- lV2 
no-' 


(2.48) 


where M is called raach number and defined as 

M = -J- (2.49) 

Solving equation (2.48) as a quadratic equation we have 

2 

tj^(t) = Fj^(t) = [t + 4 R^kd ^ )^^^ c"^ cosyj 

R **" 
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- [4 + 8 Mt (1- 

no 


R 


2 

2 


4M^(t^~4T^ + 4(R ^-h^) 

' no ' o ' 

X (l-4j\;^)"*^ 


c”^ cosy + 
c“^ cos^y)]^'^^ ^ 

(2.50) 


If we make realistic assumotion now like iu << 1 and H << R 

‘ ( 

(i*9., velocity of target is much smaller than velocity of 
Sound and depth of the channel is much smaller than range), 
then equation (2*50) simplifies to 

Fn(t) = t 4- 4R^ M c“^ cosy - 2(T^^+ 2K^ h c“^(cosY)t 

+ t^)^/^ (2.51) 

F (t) is the transmitted pulse time frame ('t ') which 
I* n 

has limits 0 ;< Fj^('t) .< T where T is the pulsev^idth. \;e 
expand ■tj^(Fj^('^)) using iviaciaurin' s series. 

1 2R k c"^cosy 

F„(t) = t-2T_^g + 4 a c-^ COSY S. t 

no 

(2.52) 

2R hi c^^cosy 

Consider t(l >• ^ — ) in equation (2.52), where 

^ no 

is the travel time of nth ray v.’hen target is stationary 

and it will be in the order of 2fv^ c"'^. Even for cosy = 1 ,. 
2R M c ’^ cosy 

. is negligable compared to '!• if M << 1. 

^no 
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Therefore, we write 


Fn(t) = t-2T^^ + 4i\ M c~^ COSY 
n no o * 


(2.53) 


where t = 0 (2R /c). 

o 

Vie now proceed to find expressions for the attenuation 
coefficients, phases and inclination angles associated with 
various rays. 


The inclination angle 0^ for the nth ray is easily 


seen from the geometry of Fig. 2.11 as 
-1 


Q =s tan 
n 


c 


P - a ^ 




J 


(2,^4) 


Spreading loss associated with nth ray is computed 
with the assumption that they on the average are cylindrical 
losses, we have for one way path from ecuation (2.28). 

L~^^^(t) = jnth ray path length] 


~ c“^lT 


2 

no 


-K H 




2R c”^ M(cosy) t 1 


-1/^ 


n' 


(2.55) 

Let be the boundary reflection loss of nth ray 

per one bottom reflection derived from equation (2.31). 

Then for N bottom reflections along nth ray we have total 
^ N 

boundary loss - [B^^] 
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Therefore total amplitude loss is given as 


N. 


\(t) = [B„] " L-^(t) I I 


(2.56) 


where is the complex reflection coefficient of the 
target. 

NoV(/ we linearize L~^(t) as follows; 


-1 


Ln* (t) = + Pv? F^^(t) + 2 Rq M F^(t) c“-^ cosy] 


“1 


-1 


« O M^F ^(t)+2R M F^(t) c“^C0SY , 

= c-2 T-^ 1 + 2 °_.-J !. ]-l 

no Y Z 

no 

(2.57) 

Again using the procedure in simplifying (2.51) we have 


( t) = (l-2fA cosy (|i - 2)) 


(2.58) 


where is the order of Fx /c, 

no o 

Hence total loss Aj^(t) is given as 


N_ 


A^(t) = [B^] " ja'^j (c T^q)"^ (l~2^a cosy(^ -2)) 


ct 

‘o 


(2.59) 


The total phase is given as 


S^(t) = [H + ^ (2.60) 

where first term gives the phase change due to surface 
reflections and second term gives the phase change due to 
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bottom reflections and€. is the phase advance due to one 
botton reflection. Here we assume the number of reflections 
and phase change 6 are not function of time as the order 
n of the ray is not large* Hence ^j^(t) is not a function 
of time, 

Knov^ing the values of and the 

received signal contribution from the nth ray can be 
computed. It is given as 

exp j | w^( t-2T^(t) ) + ], 

= ae[A^(t) exp j£%t + ] (2.61) 

where, 

«>n(t) = -2co^[T2^-Hvi2F^^(t)+2R^W F^(t)c“'^C0SY]^/2 

(2.52) 

Consider the expression given in square brackets in the 
above equation; and using the simplification procedure 
adopted in case of equation (2.51), we will have 

0^(t) = “2 Wq[T^q + M(coSY)t-2Il^ M c^^cosyJ+S^ (2.63) 

NOW we write the total acoustic field received as 
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N, 


(t) = Re[ X i\i (CT^^ )'^(1-2M coSY(|i -2) 

n o 


ct 


X exp 0) t - 2 w T «2w M(coSY)t+4u> R Cosy+S ?] 

L o o no o ’ o o ’ n- 

Rewriting the above equation as 


r(t) = Re$/[i -2M cosyC^ -2))exp j(-2w^M( cosY)t 


-1 


N 


4. 4a)^R^Mc- coSY)/x[. X [bJ ^|a^!(c 

n 

exp J(-2 Vno-^n^J * ® J j 


-2 


(2.64) 


Thus the received signal is seen to be the product of two 
terms, one of which is a function of the target velocity 
alone and the other - under the summation sign - is the 
composite signal corresponding to a target at rest. Denoting 
the later term by ®xp(j0^) we write 


N 


exp(j0 ) = E 
n 


(2.65) 


Now we have 


r(t) = [(1--2M cosy -2)) exp j(~2a)^ M(cosY)t 




« J UJ w ^ 

+ 4 w^R^ M c"-^ cosy + 0^)] e ° J 


( 2 . 66 ) 
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Comparing with the above, the r(t) given in equation (2.42) 
we observe 

=0^-2 0 )^ M cosy (t-2R^c“^) (2.67) 

and 

A(t) = Aq i-2M cosy (|^ - 2) (2.68) 

o 

\ie thus finally have, as M <<1, 

1 ^ 

r(t) 5s Re [a_ exp j(0 -2w M coSY(t-2R^c“'‘') ) ]e ° 

\j o o o 

(2.69) 

or write 

x(t) = Re [A^ exp j(+ a:^t)exp j(0Q-4w^ M cosy R^c”*^)] 

iw t 

e ° (2.70) 

where, 

(Oq = ^2 M(cosy)t (2.71) 

VVe now observe that the doppler frequency (t*^)# we 
obtained in the analysis is same as the doppler frequency 
one obtains in classical sense for a situation where the 
target is movir^ radially* 

Rewriting equation (2.70) as 

r(t) = Re [Aq exp j(0Q-4m^M cosy R^^c*^)] exp j(o)^+a^)t 


(2.72) 
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In the above equation, the term, 

Since M<<1 and T is in the order of R c”^, we 

no o ^ 

simplify it as 0^(0), 

Therefore, 

o 

I(t) = Re [S [B„] la„] {c exp(-2a>^T„„+S^)] 

exp j((o^ + u^)t (2.73) 

Therefore the composite received signal for a moving 
target situation is seen to have amplitude equal to the 
amplitude computed when the target is assumed to be at rest, 
and having a carrier frequency shifted by from trans- 
mitted frequency , where is due to a moving target. 
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2*7, Results: 

We want to investigate the variations on time 
delay differences and attenuation coefficients j a^j 
in a shallow water channel with known source and target 
depths# 

The values of the parameters we consider below are 
typical to a situation where a torpedo is homing on to 
a submarine* The following parameters are assumed: 

/ 

i) Velocity of propagation in the channel : Cj^ = 1500 m/sec* 

ii) Velocity of propagation in channel sediment: C2=1600ra/set» 

3 

iii) Density of the channel : = 1000 gm/cm • 

3 

iv) Density of the channel sediment s ~ 1250 gm/ cm • ‘ 

v) Channel depth : H = 300 m* 

vi) Source frequency : f = 30 KHz* 

vii) Transmitting pulse width : T = 100 msec, 

viii) Target velocity : V = 15 m/sec, 

ix) Direction of moving target : y = 75®, 

x) Target range : R » 500 — 1000 - 1500 — 2000m. 

xi) Target depth : Hj = 0 - 50 - IQO - 150 -- 200 - 250-3000* 

xii) Source depth : = 50 - 150 - 250. 

xiii) Beam pattern : Vertical : 40°} Horizontal : 80°. 
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xiv) Target strength : IF 70® < Y < 90® : TS 300 

else IF 0® < 10® : TS = 30 

10® < Q^< 20® : TS = 35 

20® < ©yj<45® ; TS = 60 

45® < ©^< 90° : TS = 150 

We compute the desired parameter values, for 
direct ray and the sets of reflected rays with order 
N = 0,1,2, 3, It is observed that for ^ inclination 

angle of the rays are out of beam pattern. 

The following parameters are computed for each ray 
and the equations which are used are shown in brackets, 

i) Critical angle (© ) i from equation ^z-sa.) 

ii) Inclination angle (^) J from equation C*-*^**) 

iii) Time of travel (T) ; from equation 

iv) Spreading loss (SPLS) : from equation 

v) Boundary reflection loss(BLS): from equation ^-'^9 

vi) Bottom reflection phase : from equation 

vii) Total phase introduced duo 
to boundary reflection * 

where M is the number of surface reflections, 

viii) Array factor AF : using and Beam pattern specifications 

ix) Target strength: TS: Using Yf and target strength 

specifications. 

x) Total attenuation coefficient | j : SPLS x BLS x 

AF X TS, 
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Then normalizing the parameters T and j 'a^^l with 
respect to direct ray parameters, we compute time delay 
differences and normalized attenuation coefficients { b^j 
as 



where and | 'a^j refer to the direct ray. 

The following results are observed; 

!• Fox a given range of the target and source depth, time 
delay difference of any ray path is linearly 

dependent on target depth (since the linearized version 
of the equation (2.45 ) is used). The gradient of the 
ray depends on the type of the ray. Refer graphs 2.1 
and 2 #2. 

2. For a given source and target depths, the time delay 

differences depend non-rlinearly on range of the 

target. Refer graph 2.3, 

3. Fox a given source and target depths; graph 2.4 shows 
the attenuation loss with respect to range, 

4. The spreads of time delay differences for various 

rays are shown in Table 2,1 The rays that have higher 
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order have larger spread; i.e,, if are arranged 
in increasing order, then we -observe 

^ fox *1^ < Zj 

where ^ is the spread of the delay 

This fact will jpe useful in Chapter 3, in computing 

the correlation coefficients and their probability 

densities, among signals wAiich arrive along different 

paths, 

5, The number of significant paths that contribute to the 
composite received signal at each range are shown in 
"^able 2«2* 

6, Assuming each range bin is of transmission pulsewidth 
•T^, starting from direct ray time delay 2T^, the time 
delay differences are arranged in corresponding range 
bins. This observation is helpful to see how the return 
signal energies via different paths are distributed. 

Some typical results are tabulated in Table 2 •2, 

7, Table 2,3 gives the type of rays that contribute to 

received signal at each given range. We observe that 
the number of rays contributing to the channel are same 
irrespective of source and target depths. The only 
change we see for 1000 m and 1500 m,, is due to the 
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the symmetry of ray paths with respective to source and 
target depths. This result v/ill be useful in estimating 
the arrival times as will be discussed in Chapter 3, 
f2. 1,2.2) 

8, In Graph £ we observe at a given range and source depth, 

the variations in time delay differences. The values 

taken for range (K)j source depth (H^^) and target depth 

(Hj) are arbitrary but typical, Vve observed that some 

time delay differences (t. 2 » "*5 etc.) are increasing 

with respect to target depth whereas the other time delay 

differences etc.) are decreasing with increase 

in target depth. This implies that as changes some 

paths are moving away from direct ray and some are moving 

f or 

closer. Therefore ^ a change in there is a net amount 
of change in this movement. This result will be of 
considerable importance to analyze the performance of 
conventional result in Chapter 3. 
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Table 2»2: Total number of arrivals and their 

distribution with respect to direct ray 
arrival time for different ranges. 


Range 

(R) 

Vertical 

height 

(hi 

Total 

arrivals 

Arrivals with respect 
range bin 

1 

2 

3 

4 

5 

500 m 

< 50m 

2 

2 


•» 




1. 

2 

2 

- 

- 


- 

1000 m 

< 50m 

4 

2 

2 

- 

- 

— 


> 50m 

4 

2 

2 

- 

- 

- 

1500 m 

< 50m 

6 

3 

2 

1 

- 

- 


> 50m 

6 

4 

1 


o 

- 

2000m 

0 m 

7 

■ 

3 


H 

— 


50m 

8 

H 

2 



- 


loom 

8 

H 

2 



- 


150m 

8 

B 

1 

2 

1 

- 


200m 

8 

■ 

1 

3 

mmm 

mm 


(h) : Vertical height between source/receiver 
and target depths. 
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Table 2.3: The type of rays that contribute 
at each range. 


Range 

(TO 

Vertical 

height 

(h) 

Kays that contribute 
at receiver 
(Ref, Fig. 2,6) 

500 m 

0 - 250 m 

1, 2, 3 

1000 m 

< 250 m 

2| 3, 4 


250 m 

1, 2, 3, 5 

1500 m 

< 250 m 

1, 2, 3, 4, 5, 6 


250 m 

T> 2, 3, 4^ 5, 7 

2000 m 

f 

0 - 250 m 

1, 2, 3, 4, 5, 7, 9 


(h) : Vertical height between source/receiver 
and target depths. 


CHAPTER 3 


RECEIVERS FOR TARGET DETECTION 
3.1 problem FORMULATION: 

In Chapter 2, it vjas seen that the signal received 
from a shallow water target is a weighted sum of a finite 
number of replicas of the transmitted waveform with 
different delays and a common doppler shift. The number of 
paths that contribute significant energy at the receiver are 
limited by i) the beam patterns of the receiver and the 
transmitter and ii) the number of bottom reflections a ray 
.path undergoes before it reaches the receiver. i 

In this chapter different receiver structures are ; 
considered and their relative performances studied. All the : 
receivers are aimed at the detection of a target under the ; 
usual binary hypotheses testing [l6]. ! 

i 

Let the transmitted signal be a rectangular pulse with 
carrier frequency oo^. Accordingly | 

S^(t) = f2 Re[VB^ f(t) exp(jco^t)] (3.1) ; 

I 

where i 

f(t) = C) < t ^ T I 


0 


otherwise 


(3.2) 



63 


is transmitted signal energy and the complex envelope 

f(t) has unity energy. Following our discussion in Chapter 2, 
we write the received signal component arriving via the ith 
path as 

j ( to +ax-, ) t 

Yi5(t) = {2 ae ■f(t-2T.)exp j(-2a)^T.+S..)]e ° ^ 


(3.3) 


where, 

is the attenuation suffered by the ray along 
the ith path. 


is one way travel time, 

Sj|^ is the phase introduced due to reflections at 
boundaries and target, 
and Wp is the doppler frequency. 

The composite received signal may be v;ritten as; 


Y(t) = f2 Re 


M . .. ^ j(w +oo..)t 

a^if(t-2T^)exp J(“2w^T^+S^) ]e 


(3.4) 


where M is the number of significant paths, Sj^'s are 
modelled as complex Gaussian random variables, bach of 

tu te s two parts, c ompl e x de t ermi nr sere var ra ol 
describes path losses and boundary losses, And a^^^. 


[ 

i 


3 
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zero mean complete Gaussian random variable, describes 
the reflection coefficient of the target, which we assume 
to be a slowly fluctuating point target, a^. has variance 
the phase is uniformly distributed over 

[0,27cj , 

Now we write 



The phase distribution of is seen to be uniformly 
distributed v;ith a bias component which is equal to ^jLi* ® 
phase, iVe also note la2ii ^ function of target's aspect 
and altitude angles with respect to receiver; which in turn 
depend on path inclination angles Hence a^ have in 

general unequal variances. 

'lie now consider the time delay (Tj^) and attenuation 
coefficients (a^,) normalized with respect to direct ray's 
time delay (T^) and attenuation coefficient That 


is define 
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ri A 2 (i^ - T^) 

(3.5) 

r* 

a . 

b. A 

(3.7) 


for 1=1... M. 


s arc now time delay differences and b^' s are nor- 
malized attenuation coefficients. If direct ray subscript 
is taken as '1', then 't:', = 0 and 2a3 = 1. 

J. 

Rewriting the composite received signal along with 
noise component which is considered to be a sample function 
from a zeromean complex white Gaussian random process, 
we have 

Y(t) = f2 Re[J bj_ f(t-t^)j exp^’ (w^+Wj^) tj 

+V2 Re [ 'w(t) exp(ju)^t)] (3,8) 

As we observed in Section 2*6, the coraposite return 
signal is having a common doppler shift. Therefore with 
no loss of generality we equate to zero. 

The complex envelope of the received signal is 
obtained under binary hypotheses testing as, 

Y(t) = V^t '"i ! . -i: 

W(t) ! Hq 



where 


0<t<T*’^T + Max( j) ; j 
and E[w(t) ■w(u)] = 2a „ (t ~ u) 

bT 

ue now consider various receiver configurations and 
attempt to evaluate their performances, 

3,2 OPTIMUM RECEIVER: 

The complex envelope of the received signal is given 
by equation (3.9) and is reproduced as 

Y(t) = 'b^ f(t-'t^) + w(t): 

w(t): Hq 

where attenuation coefficients ( b^) are modelled as complex 
Gaussian random variables. Time delay differences ai*e 

functions of target range, target depth and source depth. 
Hence s can be seen to be some functions of target depth 
if range and source depth are known. Since target may lie at 
any depth in the channel, may be modelled as random 
variables, m a range ^ ^ A "" ^ imin* 

Now we are interested in designing an optimum 
processor with above specifications, i.e., the problem is 
to observe Y(t) over an interval [0,T' ] and decide whether 
Hq or is true. The criterion may be either Sayes or 
Neyman-Pearson criteria. 
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Our observation is a time continuous random waveform. 
The first step is to reduce it to a set of random variables. 

l.-e observe in equation (3.9) that the received signal 
'Y(t) is a zeromean random process, with variance K^(t,u) 
given by 

K.(t,u) = E Sj, 20. E[?(t-rpf (u-ti)]+ 

'' i 

(3.10) 

'lie want to expand y('*-) using K.L. expansion method to 
obtain a set of uncorrelated random variables. 


We write 


Y(t) = 

where - 
and f 

Expanding Y('t) 


-ij, 0i(t) (3.11) 

f Y(t) 1*1 (t) 

> 

a set of complete orthonormal functions, 
by using the above functions we will have 


Yi = Sj VEt Pii + = «1 

T' . 

where f ^i 

w. = / w(t) 0j_(t) dt 

^ o 


(3.12) 

(3.13) 


and 


(3.14) 
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0i( '-) ; i = has to satisfy the following integral 

equation: 

T 

0 ^ •*■ 
where ^ ^ are eigen values* 

The above integral equation may be written asj 

T * 

E[ I qt-r.) f (u-K)]] 

J j o 

X 0j^(u) du + 2o^^ . 

Instead of solving the above difficult integral equation 
v;ith a Kernal involving the statistics of the various delays, 
' 2 ^^, we explore the following alternative approach* 

Towards that end consider the expansion of Y('t) 
terras of complete orthonorraal set of' functions ( t)^ , the 
first M of which are given by fCt-T^) for some suitably 
chosen set of values ; i = l..*M, 

V^e define, now, two vector functions as follows: 

F(t)4 ... ( 3 . 15 ) 

F(t)A [f(t-q) fCt-fj) ... f(t-rj)]'’' (3.16) 

Then we expand yCt) ovei F(t) and the function^ j^(t)J” . 

■.Ve obtain a set of random variables 
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T» 

Yi =* / V(t) f*(t~r.)dt 

o ^ 

= 'l\ bj p .(i) + 


“i = ’■‘o 

(3.17) 

T * 

where /^.(i) = / f(t-t.) )dt 

^ Q J ^ 

(3.18) 

T? 

and w. = / w(t) f*(t-'r.)dt 

0 ^ 

(3.19) 

For i = same for both the hypotheses 

and provides no additional information and hence can be 

ignored. Observe that are uncorrelated random variables. 

The cross covariance of Yj[ Yj^ given as 


(3.20) 

Nov/ we have a set of random variables Yj|^> which are 
functions parameters b. and (^.(i). he want to now formulate 

*J *J 

the problem as a composite hypotheses problem, he first 

introduce some notation. Let 

V' T 

R ^ [Vi> Y2 ••• VjviJ 

(3.21) 

B-:^ Cbj^j ^2 *•* 

(3.22) 

W 4 Cwj_, \ ... 

(3.23) 
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Accordingly we will have 

e[b] =± 0 

and EE B B"** ] =X'b (3.24) 

where A is a diagonal matrix with diagonal elenents 

equal to 2a ^ , 

^i 

Vje also define the correlation matrix, 

T' 

/ F(t) F+(t) dt (3.25) 

o 

X * 

where P.(i) = / f(t-r.) 'f*(t-r.)dt 
j 0 ^ 

The composite hypotheses problem is now to compute 


A{r ) = 


j P( R/B, f, H^)p('b ) p( PJ d Bd ‘ 


P (R /H;, ) 


(3.26) 


Towards evaluating A (R ), we first consider evaluating 

A( ii/p) by averaging A (H /B , P) over the space of B 

parameters , A (R /P) is a complex Gaussian vector under both 

Hs and H, and hence its solution is obtained by computing 
U 1 

the conditional likelihood ratio, which is given in 


general as 
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A ( Vp) = 


expC- I ( Qj_( R - Mj_) 

exp [• ri^ - M^) (R - M^)] 


where, 

E[ R/Hj^] = E[ R/HqJ = 1% 

E[( R-M^)(R^- M^'^)/Hj^] = Kjl 


(3.27) 


E[(R - M^) (R+ - Mo'^)/Ho] = 

'h = V" » = ^o"^- 

In the present problem, we have 


R ^ ■p'^ B +li i 


S9 

W ; 

Ho 

(3.28) 

a nd 

o 

fl 

o 

!1 




= £t b 

P “t 

(3.29) 

^0 

= IS/ = 2a/ 

I 

(3.30) 


Then the conditipnal likelihood ratio is simplified to 


/\ ( R/P) = 1 R (Qq “ % ^ ^ ^ 

^1 


(3.31) 
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Now we need to compute to obtain a simple 

structure for the above ratio. 

Consider the signal component of covariance matrix 
Kj^, and it is given as 

^hS^^t^^^b ^ (3.32) 

This is a non-diagonal matrix implying that the compo- 
nents of vector K are hot independent, Vve now diagonalize 
the matrix by solving its characteristic equation and 

obtain the eigen values and eigen vectors. If is not 
positive definite, we augment the eigen vectors with some 
more to form a complete set. V/e have seen to be a 
diagonal matrix with equal diagonal elements (20^^ ) and 
so the diagonalizing matrix (^) will not effect the matrix 
K So write 

'% = ‘hs 

and (3,33) 

If 2cf ^ are the diagonal elements or eigen values of 
they represent the energy components of received signal 
along different basis functions and also their sum is 
always constant irrespective of the choice of 0 . Now 

» - Q ' is obtained as 
JL* O 
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Substituting the above matrix in conditional likelihood 
ratio, we obtain, 


(3.34) 


A( R'/P) = exp (3.35) 

i “i +”10 2 “a. 




2.^2 ^i 


a. 

10 ) 


where = jj**" R 


(3.36) 


The above ratio is developed with apriori knowledge 
of P i.e,, P.(i) for i| l.,.M. Now if we have the 
probabilistic knowledge of P, then the optimum processor 
would be 

(r' ) = / ( r'/p) p(p) dp 

p 

This integration in general will be difficult and we 
consider a case where M = 2. 
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where A a E^( ^ ^^2 p^(^i)) 

K22sA® (2) + 0^1 /^l^'^') • 

and ^215 “ •'‘ 12 &A *'— P j_i.^) ( x(2) + 

Solving the characteristic equation we will have 

2 r. 

1 * A“^ and ;i 2 = ® 

if c^ < (A-B)^/4. 
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then the conditional likelihood ratio is given as 


A(R’ /p) = 



1 2 2a/ ^ ^*a/ 





and then 


/5(r') * /A ( rVp) p(P) dP^ (3.37a) 

where P. (i) * / f(t-7". ) 'f*(t-T. )dt. 

J o ^ ^ 

Now pi(i) is a random variable and it's probability 
density depends on the density of Tr.. From • .(i) we now 

%J *J 

have to compute joint densities and 7)2 to evaluate the 
integral (3.374. 


This approach of evaluating likelihood ratio involves 
extensive computation and the receiver implementation would 
also be difficult} even for the case M = 2. 

But however, if we could compute the joint densities 
of and ^2 assume high signal to noise ratios, we may 

be able to obtain a simpler structure. The eigen values 
and p\ ^ are the received signal energy components and if 
> i >> in equation (3,35), we see that 


A ( R*) = 2a exp 3 ( )] ^ t) (3.38) 

w 2d 2 

w 
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where 


a 





( 3 . 39 ) 


Likelihood ratio is obtained as. 





* 


The performance of this receiver can be easily obtained as 
^ (R* ) is having central chi-square distribution, 

l/^e can also observe that for high order M, the 
sufficient statistic is still same i,e., 

^(R') = S“=l ^ Y* 

i 

if high SNR condition is satisfied. But we need to compute 
Mth order joint density function and then evaluate the 
corresponding integral. 

For low signal to noise ratios, the receiver may not 
be as simple as shown above. By appropriate modelling of 
densities ofX^t we may still be able to realize the optimum 
processor, 

We now consider some other receivers wherein we have 
precise information of the time delay parameters. 
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3,3 OPTIivUM RECEIVER VVHEN TIME DELAYS ARE PERFECTLY KNOV^: 

In Section 3*2 we have seen an optimum receiver 
foxftulation when B and p are random vectors and matrices 
respectively. This was found to be a difficult problem even 
for M=2* Now we consider an optimum processor for the 
situation in which time delays ' l- ^>^^1 therefore p are 
perfectly kno\A/n* The performance of this receiver if 
evaluated gives a bound on the performance of the receiver 
discussed in Section 3,2* In Section 3*4 when we discuss 
a conventional receiver with perfectly knovm delays, the 
performance of the optimum receiver will give us an estimate 
of the degradation that will occur w;ith a conventional 
receiver* 


We as seen in Section 3,2, again consider a set of J 
orthonormal functions^T^^( t) ^ wf which the first M functions; 
will be F(t) (eqn* (3*15), This choice of F(t) is essential { 
to obtain uncorrelated noise components when yCt) is expanded j 
Under known P, zeromean Gaussian under both the j 
hypotheses, but will be non-diagonal. Solution of it’s j 
characteristic equation gives a transformation matrix which i 
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This is to be compared against a threshold y* where 

2 2 

yif- =^lnT) + In —4.. ^ ^ ^3 

a 

0 ) 

The receiver may be implemented as shown in Fig, 3.1. 
3,3.1 Performance Evaluation: 


In the hypotheses testing problems we know that the 
log-likelihood ratio is given by 1; ( R) and is compared 
against a threshold for detection of target. R) is a 
random variable and it's density depends on the hypotheses 
considered. If the densities of ^ ( R) are known under 
and in a binary hypotheses problemj we can compute the 
false alarm and probability of miss obtained for a given 
threshold and signal to noise energy. 

As seen in eqn, (3,^o), the random variable £( R) is a 
sum of squares M Gaussian random variables with unequal 
weights. 


The computation of py^^ ( L/W^) is a difficult job 
as it involves unequal v/eights for the variables. Therefore 
we compute the performance bounds for the above receiver 
through moment generating functions [l7]. These bounds 
are called chernoff's bounds. 


These bounds are found by computing a function 
^(s) given as, 
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}i(s) 


oo 


(3.42) 


where 's* is taken as a parameter such that s > o. Then 
false alarm is given as [l7], 

,2 


Pp < Exp|j[p,/s)-.s |a(s) -f p.{s)) erf c*(syp,(s) ) 


^ 'S > 0 


:3.43) 


and probability of miss is given by 


2 

exp(ji(s)-(l-s) p(s) + ^(s)) erfc*((l-s)y ii(s)) 

s < 1 (3.44) 

where p.(s) and ^(s) are the first and second derivatives of 


Vi(s) 


The minimum bound is obtained by computing 


given as, 


^l(s)l 


= T) 


s»s 


min 


(3.4b) 


where t) is the threshold. For equal probability of occurrance 
and equal cost functions, we have t) = 0. 

Computing a satisfying the above equation, now 

false alarm and probability of miss are given asj 
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Pp < Exp [p,(s) + |~ ^i{s)]erl c*(sVtr(s) ) j .'-z< 4 ^) 


s=s . ^ 

min 


2 

PjVi S. [p-ls) 3 erfc*( sf^‘(s) ) j 

where 0 < s 1, 


s=s . 

ram 




erfc’> (s . V" Vi(Sn,in^^ — 

min min ^f2% s . fu(s . ) 

^ min'^' min^ 


Similarly if (l-s . ) fiiis . ) > 3 then, 
" ' min' min * 


-s^. /2 

min'" 


erfc* ( (l-s)Yp(s)) I 


exp( C . V:| - } , ■'il(s)) 




Therefore, 

p < — 1 exp [^i(s)3 j 

f2%s^ *lICs) ®~®min 

p < exp t^i(s)] 1 (3.49) 

■f2ix(l-s)‘^ p(s) min 

3,3.2 To Compute Performance Bounds for Equation (3.4-o ): 

’We have seen in Section 3.3.1that p(s) is given as 
€0 

p(s) = in / [A(R')]® P„,/H ( R/H^)®'. 

..00 ' ' O 
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We know A (R«) and ( R/H^) where A (R* ) 

.s given by egn* (3,38 ) where is given as 

p” 2 

i ^1 n 

0 


Qi - Ci' = 




CT, 


IL 


/ 2 2x 2 

( a, +a )o 

'■ M 0) ' OJ 


Then ^l(s) can be evaluated and is seen equal to 


M 


p(s) = In / [ 7c ( 


^—2)3 

i=si cr, +cf 


% ,^s/2 


■s 


<> n nrw 
1 0 ) 


Ti [exp [| + 2 

i=l o. 


w 




.2 


" IT'TS 7:~~2 -'72 J (3.50) 


“i 


f2tta 


n 


2 a. 


(j) 


Completing the squares in exponent and solving the integral 
p(s) is simplified as, 

, („ 2+<, 2)l-s 

t.{s) = L i in [ ; " ■■■ ] 

(l-s)a. + 


Now to compute Chernoff’s bounds, we proceed as 
follows: 

i) Equation (3*5© ) is computed for various values of 's‘ 
in the interval [0,l]. 
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ii) The minimum value of ij.(s) is found and corresponding 
’s' is the and at this 's' we have p(s) = 0. 

iii) IVe now compute p(s) as 


p(s) 


1 , r '“l /‘’u. ) 

2 i^i L 77;“,— ;“ - 2 


(l+(l-s) 


iv) Calculate u(s . ), 

min-^ 

V) Define x A and 

now if erfc* (x) > 3, the false alarm is given by 
eqn. (3,48 ) or else by eqn. (3 ,a 6. ), Similarly if 
erfc*(y) > 3> the probability of miss is given by 
eqn. (S.-o.-'! ) or else by equation (3. -*. 7 ). 

vi) The total probability of error is given by, under equal 
likely conditions, Pj,(e) = 2 % 2 


The results on these bounds are discussed in Section 3,5, 
i/Ve now proceed to consider a much simpler processer called 
conventional receiver. 


3,3,3 Selection of Transmitted V^aveform to -Reduce Receiver 
Complexity; 

The complexity of the receiver structure considered 
in Section 3.3 arises from the fact that the signal contri- 
butions reaching the receiver via the various ray paths, are 
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not orthogonal. One way of making them orthogonal would be 
to choose a transmitted wave form f(t), such that its 
correlation function 

T 

Rf(T)A / f(t) f"(t+r)dt (3.51) 

0 

satisfies the property, 

Rf(r) = 0 for all f > 

where, 

T < Min -r.-l ) (3.52) 

It may be recalled that all i^j, are non-negative 

differences in delays, between the reflected ray paths 
and the direct ray, 

A linear frequency modulated waveform or a Barker 
coded sequence of pulses or a PR3S sequence, with appropriate 
choice of parameters will approximately satisfy the condition 
given in equation (3,52 )• 

Assuming as in the previous section that the time delays 
T'- are perfectly known at the receiver, with a transmitted 
waveform f(t) chosen as above, we can easily show the simpli- 
fication that will occur in the receiver structure, compared 
to the receiver shown in Section 3,3. For this purpose 



84 


instead of conside.clng the subset of the orthonormal functions 
^ 2 .)* i=l«.»ivl as the first M members of the COM set 
Section 3,2, we nov^i choose them to be 

/ ■' 1=1 

'i/ith Such a choice in view of equation 
(3,52 )» it is easily seen that 

Pj(i) = ^ 

0 otherwise (3,53) 

As before the hypotheses testing problem simplifies to 


= Zj 5 "i ' M 


w. 

1 


H. 


(3.54) 


The covariance matrix of the measurements is now a diagonal 
matrix with i,jth element given as 


eCy 


1 '3 


] = 


2E, 


0? + 2cf ^ 

D. (!) 

1 



Proceeding as in Section 3.2, the likelihood ratio is given 
as 


iVJ* 

A ( R) = IX [ 


^ j exp [ 2. ^ 

i=l E^o;: +aj ^ 2a/(£^o ^+a 


E a ^ 

t%i 11 2n 

2, ””“2 N 


't“b^-''w 


w 


The log-likelihood ratio isj 


(3.55) 


£( R) 


^t 




Yii 


(3.56) 
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and the threshold y 


[inn ~ In 


M 

i=l 





to 



The receiver implementation is shown in Fig. 3.2. 


The performance bounds of this receiver can be computed, 
his can be easily carried out with a substitution of 


t b 


/(^t^^b^ + in place of (3«50 ) 


t"b. ‘ ''w 
Therefore ij,(s) is given as 




2 xl-s 


^(s) 


Z In 


(l-s) £. o, ^ + o ^ 
' t b^ w 


and 


E. cr._^ /a, 


|a(s) 


. 't ^b, ' ''0) o 

1 2 [ — 5 — -“ 5 - r 

l+(l-s)CE^aj^2 ^^2 


We can compui:e Chernoff bounds by computing 
and then substituting this in corresponding equations. 

The results obtained by this method are discussed 
in Section 3.5, 

3.4 CONVENTIONAL RECEIVER: 

Consider the equation (3,9} in which we assumed b^ are 
zero mean complex Gaussian random variables , s, the time 
delay differences are assumed to be known. Depending on the 
values of T. andr.^ we have f(t-%) and f(t-Xt) overlapping. 

1 J ^ 
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The conventional receiver is designed v/ith the 
/iew that only o^rec'c ray and the white noise component 
are present, ihen knowing the direct ray time delay (Tj^), 
we compute 


Yi = / y( f ("t - T, ) .dt 

a- JL 


( 3 . 57 ) 


or equivalently when we assume = 0, we have 



yCt) 


■5^ 

f (t)dt 


Ei b. /0^(1)+ ; Hj_ 


^1 * ^0 


( 3 , 58 ) 


vjhere, 

i(l)= y f (t- i) ftt- ^)dt (3.55) 

and-,(i) .are knovm. 

X 

Therefore is a sum of utmost independent zero mean 
complex Gaussian random variables along with a white noise 
sa-i;ple, under Hence Yj_ is a complex Gaussian random 

variable under with zero mean and variance 2{^i ) where 

2 

Under il„, y, is having variance 2a . Then 

OX ^ 

likelihood ratio is; 
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A ( '\) = 2 [(Yi)' 




”2 2 ' 2 ^ ^ ' 
a, (0, +0 ^ 

CO '' 1 (,0 


and sufficient statistics is 


n R) = IyxI^ ^ Y* 


Where, 


= [in'o - In 


(0 _ 

« 2 
°1 


„ 2 ,„ 2 ^„ 2, 

0 ( CT, +0 ) 

CO ' 1 CO ' 


Perfor.nance of conventional receiver ; 


We have 


•£ ( H) = |k,|2 = If ^(t) ^)dt|2 ( 


and we know the receiver performance is given by, 


A 


E[(Rj^)^/% -. E[( Rj^)^/Hp] 

E[1 RiI^/Hq ] 

^ ^ 1 (1) 

• ^3- from equation (3*60) 


CO 


Mow vae know that false alarm Pp and probability of 
Pq are related asj 

Pp = 


3,61) 


3,62) 


3.63) 


3.64) 


3.65) 


detection 


(3.66) 
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Now for various values of SNR., the performance of conven- 
tional receiver can be computed. In Section 3*5 we compare 
the performance of this receiver with the performance bounds 
obtained for the optimum receiver with known time delay 
differences considered in section 3,3. 

3.5 PSRFOhiviANCE CALCUL^^TIONS OF CONVENTIONAL RECEIVER s 

The performance of a conventional receiver can be 
evaluated from the plots of probability of detection (P^) 
with respect to SNR. for given values of probability of false 
alarm (Pp) . The relation between Pp and P^ is given as 
(from equation (3.66); 

Pf=(PdV'^'‘‘ (3.66) 


where. 






9 

Here a. is the variance of the reflection coefficient which 
^i 

also includes path and boundary losses, Pj^(l) is given as 
T 

PAD = / f(t-'r.) f*(t-'£'.)dt 

o 

p 

he observe that 2E^ to be the total energy that arrives 

via ith path and therefore 2E^ P^^(l) is the amount of 
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energy that ith path contrioutes in the performance factor 

'A’* signal to noise ratio, then 

2 2 

cfj^ /^7(l) gives the multiplicating factor, which is the 
i 

improvement in the amount of energy 

that arrives in range bin 1, This multiplicating factor, we 
call as multipath effect (Ivl.h.). Therefore, 


aI « ij. # 




1 

As we observe from equation (3.65 ), the magnitude of M.E. 
-has direct bearing on i-iOC curves, computed for a receiver 
which receives only direct path. As M.E. > 1 always in a 
multipath environment, the receiver perforins better than the 
receiver which is operating in single path environment. 

Hence in the following v/e explore the characteristics of 
channel which effect M.E. and the performance of the receiver 


I 


Define 'h' to be the vertical cdstance between source i 

i 

and target i.e., hA where is source depth and Hj I 

is the target depth. | 

L 

The magnitude of M.E, varies with respect to range, | 

source and target depths, we plotted various graphs ! 

I 

depicting the variations of M.E. with the above parameters, i 
I've infer the following facts. I 
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This graph is plotted for the variations in M.E, 
with respect to range for different h, with Hp, = 150 ra. If 
the source is located at the middle of the channel, the Jvi.b. 
linearly increases with range irrespective of the target 
depth. But for a given range, the M.E. increases with h. At 
h = 0, the M.E. is minimum. This is due to less number of 
arrivals’ into first range bin (BBl) compared to h;>0 (refer 
Table 3.1 ), 

we observe that the number of arrivals in RBI decides 
the initial value of M.h. and iu.E, increases with range and 
this is because the net amount of overlap in RBI by other 
rays, is increasing. 

Gr aphs 3.2 to 3.3 : This graph indicates the changes in M.E, 
with range, considering ’h' as a parameter, with = 50, 250m. 
If h is a maximum, i.e., in a situation where source is close 
to one boundary and the target is close to another boundary, ' 

the M.E. is negligible at short ranges and has steep gradient i 

! 

With respect to range, he also observe, with a decrease in : 

i 

h the gradient also is reduced. For h equal to 200m, the j 

i 

increase in M.E, is 20Q/i over a range of ISOOra, and for | 

h equal to zero, the increase is just 10/- for the corresponding- 

! 

j 

range. But it is seen that, for h=0 at short ranges, the i 

M.E. is 70/» higher than the M.E. at h=200m. These facts can be ! 
explained as follows, i 
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For h = 2CDm, in Table 3,1 we observe that three 
new paths arrive in liBl when target range changes from 500 
to 1000m, This contributes to the steep gradient corres- 
ponding to j this curve. Similar effect, but to a lesser 
extent is seen for other 'h' and hence their gradients are 
less. For h=200mj there are no new arrivals after ICOOra, 
but because the net amount of overlap increases in liBl with 
respect to range, the Ivl.E. also increases. For h = 0 m., the 
net amount of overlap is not changing raarkedly and hence 
there is no relative improvement beyond 100m range. At 
range of 500m, the increases inversely proportional to h, 

because the amount of overlap of the two rays present in RBI 
increases with respect to h. (Refer Section 2,6), 

Gra ph 3,4, 3,5 : This graph indicates the variations in M.L, 
with respect to target depth for a given Hp^, taking range 
as varying parameters for different curves. 

When the source is at the middle of the channel, then 
M.E, is almost constant with target depth irrespective 
range. These curves only show a small reduction in iVi.B. when 
h=0. This is due to number of arrivals at h=0 are less than 
number of arrivals for h > 0 (refer Table 3,1 ). If the 
Source is not at the middle of the channel, M.C. increases 
with h. For range 1000, the k.E. initially decreases with 

t 
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h, because there are no new arrivals and the only one 
is present and this increases with h (refer Table 3.1 ). For 
range 500m, the ivi.b, decreases inversely proportional to h, 
because there are only two rays in RBI and as ‘h' increases 
the second ray moves out of RBI. 

Now we propose the follov/ing conclusions; 

ivi.b. is a function of both target range and *h' ; in 
other Words ivi.E, is a function of number arrivals in RBI and 
the amount of net overlap in KBl. .ie also observed from 
equation ( 3,66) that is directly proportional to k.E,, for 
a given and SNR, Vie now propose that if the source is in 
the middle of the channel, the il.E, is minimum and therefore 
the improvement on Pj^ is also proportionally less. If the 
Source is not in the laiddle of the channel, the detection 
probability is markedly larger at long ranges as h increases 
and reverse the case for short ranges. Fox a given source 
power, Si'iR decreases with range and we can effect a reduction 
in Source power to an order of 6 dB if the target is detected 
with source located at maximum possible h at long ranges. 

For short ranges irrespective of raxoge, we have approxi- 
mately 2 dB reduction due to M.E, 
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3,5,1 Comparison of Conventional receiver performance with 
performance bounds of optimum receivers in Sec. 

3.3 and 3,3.3; 

We computed the performance bounds of tho receivers 
Suggested in Section 3,3 and 3*3,3 i,e, for receivers where 
delay times are perfectly known. The receiver in Sec. 3,3 
is having rectangular transmitting pulse, where as the one 
in Sec, 3.3,1 is having modulated transmitting pulse which 
will result in mutually orthogonal signal waveforms along 
different paths. The performance bounds are calculated 
following the development in Sec. 3,3.2. For the same false 
alarm the probability of miss for conventional receiver in 
Sec. 3,4 is also computed. The results are tabulated in 
Table 3,2, The parameter values used are arbitrary but similar 
results hold good for all sets of parameters of interest. 

The following inferences are made: 

1, The degradation in conventional receiver cannot be exactly 
found as the performance results corresponding to the 
optimum receiver are bounds. For a given Pp and SNR, the 
of conventional receiver is always less than that of 
optimum receivers. Hence the conventional receiver in 
Sec, 3.4 offers satisfactory performance. 

2. The conventional receiver designed for the receiver 

discussed in Sec. 3,3,3 Vi/ill offer no improvement over the 
conventional receiver designed then the channel is thought 
to be unbounded. 
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CHAPTER 4 


CONCLUSIONS 


' In this thesis we made an attempt to model the shallow 
water targets, and then we investigated some receiver 
s cructures for detecting them under white noise conditions 
by an active SOD.-\R. 

Towards this, we have used ray theory to develops 
the acoustic transmission in shallow waters and it resulted 
in a discrete multipath structure. Each ray path is 
characterized by parameters 5 travel time, spreading loss, 
and boundary loss and phase shift. Hence the received 
signal energy is the superposition of signal energies 
arriving along different paths after getting reflected from 
the target, he have seen the number of ray paths that 
contrioute significant energy at the receiver are limited by 
SODAR beam pattern and critical angle associated with bottom 
reflection, he then analysed the variations in total loss 
suffered by a ray and it‘s travel time with respective 
different ranges of target from receiver, for source depth 
and target depth. 

For an iso-speed channel, assuming target velocity 
to be much smaller than sound velocity and channel depth 
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to be much smaller than direct ray range, we obtained the 
linearized model for the received signal for a SODAR from 
a moving target. The received signal energy was interpre- 
ted to be same as that of a stationary target with doppler 
shift in carrier frequency 'independent of number of rays 
that contribute at receiver. This common doppler shift 
for the composite signal allows us to look for a simple 
receiver configurations. 

The detection problem for the shallow water target by 
bODAR is formulated, i/iie assumed the combined effect of 
path losses, boundary losses and losses due to reflection at 
target as attenuation coefficient of the signal arriving 
along the path and it is modelled as zero mean complex 
Gaussian random variable. If the arrival times are also 
treated as random variables, then the optimum receiver for 
this characterization is practical only at high signal to 
noise ratios. Assuming perfect knowledge of the time of 
arrivals, the optimum, receiver v-jas investigated and analysed. 
The conventional receiver for this characterization found 
to, be adequate. Usage of special signal waveforms (Barker 
coded, LFM etc.) for transmitter pulse resulted in simpler 
optimum receiver structure. But the conventional receiver for 
this specifications gives no additional advantages. 
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Suggestions for Future Work: 

As an extension of the work done in this thesis an 
appropriate model for the channel where the sound velocity is 
a function of depth and the boundaries are of complex 
structure. The target may also move along any specified path. 

The model may also include the effects of reveroaration 
which then will be more realistic. An appropriate model 
for the arrival to times may be attempted to obtain an 
optimum receiver in a composite hypotheses situation. 

For the optimum receiver analysed in Section 3,2, wo 
assumed to have perfect knowledge of time delay differences. 
Some efficient estimation methods may be investigated. One 
possible way of estimating travel times is as so follows: 

With an apriori knowledge of approximate range and 
receiver depth, the time delay differences (r^) may be 
modelled as uniformly distributed random variables. If 

T . If T is the correlation time of the 

' 1 max 1 min. o 

transmitting signal and KT^ , then using a scheme shown in 

Fig. 4.1 we may estimate time delays. 
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